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Abstract: In species of similar shape and size, biometric analyses make it possible to establish dif- 
ferences. Within one species, biometrics can help researchers to detect differences between popula- 
tions and analyze their adaptations to environmental conditions. Until now little was known about 
the biometrics of the Iberian populations of Boyeria irene (Odonata: Aeshnidae), a large species 
living mainly in southern Europe. Eight biometric variables were studied in male and female exuviae 
of B. irene collected in three rivers of the Iberian Peninsula, with the objective of ascertaining which 
are the most suitable populations of this species to study. An analysis of principal components (PCA) 
shows that lengths of the epiproct, paraproct and prementum are the most influential in each of the 
three populations. The other variables (head width, body length, length of the gonapophyses, maxi- 
mum and minimum width of the prementum) proved not to be relevant in this context. 
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Introduction 


Biometric analysis makes it possible to describe interspecific differences that facilitate the separation 
of larvae of similar insects on the basis of their appearance (Landwer & Sites, 2006), and allows us 
to analyze the adaptations of some populations to environmental conditions (Hassall et al., 2008). At 
intraspecific level, a biometric analysis makes it possible to observe differences that would otherwise 
pass unnoticed. Thus, for example, the size of certain body structures in larvae of Odonata varies 
significantly between populations of Anax imperator Leach, 1815 living geographically very close to 
each other (Casanueva et al., 2017), and between populations of Cordulegaster boltonii (Donovan, 
1807) living far away from each other (Casanueva et al., 2020). Whether these characteristics coin- 
cide or not with different genetic traits has not yet been proven in many species, since there is coin- 
cidence in some (Pilgrim et al., 2002) but not in others (Huang et al., 2020). Today, modern morpho- 
metric techniques allow the analysis of body structures with more precision. Establishing whether 
and how these structures vary (Klingenberg & Marugan-Lobon, 2013) will help us to better under- 
stand evolutionary processes. 


The larvae that emerge from the water leave behind the corresponding exuviae on the riverbanks 
where rheophilous odonates live (Corbet, 2004). The study of these exuviae presents numerous ad- 
vantages as they are easily obtained in large numbers without excessive effort, their species can eas- 
ily be identified (Doucet, 2010) and they remain in the same place for a long time (Aliberti Lubertazzi 
& Ginsberg, 2009). 


There are numerous variables that can be measured in exuviae and larvae (Doucet, 2010), and their 
selection will depend on the species being analyzed. So, for example, the lateral or back spines are 
measured for the purposes of studying their antipredator function (Sahlén et al., 2008; McCauley et 


= 


Escola Superior Agraria de Castelo Branco (ESA-IPCB/CEABN-INBIO), Qta. Sra. Mércules, 6000 Cas- 
telo Branco, Portugal. 

Departamento de Ciencias Experimentales, Universidad Europea Miguel de Cervantes, 47012 Vallado- 
lid, Spain. 

Universidad Católica de Ávila, Calle Canteros s/n, 05005 Ávila, Spain. 

Departamento de Biología Ambiental, Facultad de Ciencias, Universidad de Navarra, 31080 Pamplona, 
Spain 

* corresponding author: Email: pcasanuevafduemc.es 


N 


Rw 


8 OpenAccess. © 2021 Luisa Nunes, published online 7 Jul 2021 by Wachholtz (O89 ° 
This work is licensed under the Creative Commons Attribution 4.0 Public License. 158 


@ Wachholtz 


Nunes et al.: Useful biometric variables in Iberian exuviae of Boyeria irene 


al., 2008; Leipelt et al. 2010; Johansson et al., 2017), while their prementum length can be measured 
to study their effectiveness in capturing prey (Lagrue et al., 2011). 


Boyeria irene (Fonscolombe, 1838) (Odonata: Aeshnidae) is a large species (total length of larvae 
3-4 cm in the final stage of development, designated as F-0), which makes finding its exuviae easy. 
Itis distributed mainly 1n Mediterranean Europe (Boudot et al., 2009), where it occupies the rivers of 
mountain ranges, often associated with a few species of other anisopteres such as Cordulegaster 
boltonii, Onychogomphus forcipatus (Linnaeus, 1758) and O. uncatus (Charpentier, 1840) (Ocharan 
& Torralba Burrial, 2004). It is a semivoltine or partivoltine species (Ferreras-Romero, 1997; Velasco 
et al., 2018), with a pattern of emergence suitable for a summer species (sensu Corbet, 1954, 1964). 
Its larval cycle was analyzed by Ferreras-Romero (1997), who established several stages of develop- 
ment combining the head width with the length of the wing sheaths. The morphology and biometry 
of B. irene larvae in Europe has been analyzed by Wenger (1959), Wildermuth (2005) and Müller et 
al. (2012), among others. Müller et al. (2012) proposed five measures (lengths of body, abdomen, 
cerci, paraproct and epiproct) that distinguish the larvae of B. irene and B. cretensis. In the Iberian 
Peninsula very little is known about its larval biometry, except for the work cited by Ferreras-Rome- 
ro (1997) related to larval stages, and the length of exuviae provided by Prunier (2005). Nothing is 
known about the possible geographic variation with respect to its biometrics in Iberian populations. 


The aim of this work is to analyze some biometric variables in B. irene, usually taken in exuviae, 
with the objective of ascertaining which are the most appropriate ones to study Iberian populations. 
This is a preliminary step to the analysis of the biometric variation of these populations, which can 
show differences and similarities, as 1s the case with Cordulegaster boltonii (Casanueva et al., 2020). 
In our case, we established two conditions that the populations chosen for study must meet: 1) they 
must be geographically isolated; 2) the values of the variables analyzed must differ between popula- 
tions, to ensure in this way that a wide spectrum of sizes is covered. For the variables to be useful in 
a later study on a wider geographic scale, there should be two or three variables that are representa- 
tive within each population, taking into account gender-related differences in size. Furthermore, in 
order for these variables to be used, they must be the same in the populations analyzed. 


Material and methods 


Exuviae samples were gathered in three rivers of the Iberian Peninsula, chosen for being in zones of 
different climatology and altitude: one in Portugal (River Cova) and the others in Spain (Rivers Frío 
and Eresma). Their characteristics are summarized in Table 1. The geographical distances are 171 
km between Cova and Frío, 228 km between Frío and Eresma, and 325 km between Cova and Eres- 
ma. Sampling was carried out in June, 2013 (Frío and Eresma) and in June, 2015 (Cova), by carefully 
inspecting the substrates (ground, emergent vegetation, stones, tree trunks, logs, etc.) located in the 
river banks. 


Table 1. Some characteristics of the three sampled rivers. Tm: annual mean temperature. R: annual 
precipitation. Tm and R from Ninyerola et al. (2005). 
































Cova Frío Eresma 
Municipality Ribeira de Pena Villasrubias San Ildefonso 
Coordinates 41.5821°N 40.3154"N 40.8488"N 
7.7847°W 6.6347°W 4.0179°W 
Altitude (m a.s.1.) 495 830 1200 
Tm (°C) 12.9 12.6 10.4 
R (mm) 1009 1044 835 
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The following variables were measured in each exuvia according to the criteria established by 
Wenger (1959): head width (HW); body length (BL), obtained only on straight exuviae, without cur- 
vatures produced during ecdysis; paraproct length (PA); epiproct length (EP); gonapophyses length 
(GO, only in females); maximum prementum width (W+); minimum prementum width (W-); pre- 
mentum length (PR). They are all expressed in mm and were obtained by the same person (FC) with 
a stereo microscope equipped with a millimetre eyepiece, except for BL, measured with a 0.01 mm 
precision digital calliper. 


No attempt was made to measure the length of the cerci because B. irene presents two morphs: 
brachycerca (short cerci) and typica (long cerci) (Wenger 1959), which can only be observed in fe- 
males and this dimorphism varies widely between populations (Mikolajewski et al., 2000; Ferreras- 
Romero & Suhling, 2021). For each biometric variable, normality tests (Shapiro-Wilk) and homoge- 
neity (Barlett’s) test were performed. When these conditions were fulfilled, variables were compared 
based on genders and among rivers through ANOVA, and when they were not fulfilled, they were 
compared using the Kruskal-Wallis test. These calculations were made by STATA 12.1 software. 
Furthermore, an analysis of principal components (PCA) was carried out to establish which variables 
identify each population best, and those variables were logarithmically transformed. Likewise, the 
three studied populations were compared by means of another PCA. These calculations were carried 
out using the Past 3.15 pro software (Hammer et al., 2001). 


Results 


A total of 356 exuviae were collected, 218 in River Cova (94 males, 124 females), 83 in River Frio (30 
males, 53 females) and 55 in River Eresma (21 males, 34 females). The average values and standard 
deviation of measured variables are included in Table 2. 


Table 2. Mean values t SD of the variables obtained in male and female exuviae of Boyeria irene from the 
Cova, Frío and Eresma rivers. In parentheses, range of values. HW: head width. BL: body length. PA: length 
of the paraproct. EP: length of epiproct. GO: length of the gonapophyses. W+: maximum width of the pre- 
mentum. W-: minimum width of the prementum. PR: length of the prementum. All measurements in mm. 

















































































































Male Male Male Female Female Female 
Cova Frío Eresma Cova Frío Eresma 
HW |8.07 € 0.16 7.79 + 0.20 7.954 0.15 8.33 +0.16 8.08 + 0.31 8.21+0.15 
(7.65 — 8.40) 1(7.80—8.25) |(7.70—8.20) 1(7.80— 8.70) | (7.30 — 9.00) (7.90 — 8.55) 
BL |36.24+1.15 36.75 + 1.15 37.80 + 1.28 37.91 € 1.14 37.08 + 1.55 39.15 x 1.15 
(33.10 — 38.94) | (34.26 — 39.12) | (35.41 — 39.88) | (35.08 — 41.07) | (34.17 — 40.51) |(37.41 — 41.94) 
PA 13.92-40.14 3.822 0.15 3.83 + 0.13 4.00 + 0.13 3.86 + 0.25 3.93 + 0.18 
(3.60—4.40) |(3.50-4.15) |(3.60—4.05) 1(3.60—4.25) |(3.30 — 4.40) (3.60 — 4.30) 
EP |2.88 +0.16 2.70 + 0.19 2.92 + 0.13 2.84 + 0.14 2.71 + 0.24 2.85 +0.14 
(2.50 — 3.35) |(2.30-3.10) (2.70-3.25) |(2.50-3.15) |(2.20- 3.50) (2.60 — 3.20) 
GO 3.36 € 0.11 3.33 + 0.17 3.51 +0.11 
(3.00 — 3.60) | (2.90 — 3.60) (3.25 — 3.80) 
W+ (4.28 € 0.11 4.11 +0.13 4.16 + 0.13 4.43 € 0.10 4.21 € 0.17 4.31 € 0.12 
(4.00—4.65) |(3.90-4.35) |(3.90— 4.40) |(4.20-4.70)  1(3.90— 4.50) (4.05 — 4.55) 
W- [2.33 2 0.07 2.25 + 0.09 2.29 + 0.07 2.41 + 0.08 2.27+0.12 2.38 + 0.06 
(2.15—2.55) |(2.00-—2.45) |(2.20-—2.40) |(2.20-2.60)  1(2.00— 2.50) (2.30 — 2.55) 
PR |648+0.21 5.94 + 0.43 6.22 + 0.19 6.69 £0.18 6.14 + 0.42 6.43 + 0.17 
(5.90—6.90) |(5.12-—6.70) 1(5.70—6.60) 1(6.05— 7.10) | ((5.40— 7.25) (6.15 — 6.80) 
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All variables were adjusted to normality and homogeneity, and all were significantly larger in females 
than in males, except for the epiproct in the River Cova (Table 3). In the River Frio, all variables were 
adjusted to normality and homogeneity except HW, BL and PA. In this case, variables HW, W+ and 
PR were significantly larger in females than males (Table 3). In the River Eresma, all variables were 
adjusted to normality and homogeneity, except BL. All variables were significantly larger in females 
than in males, except the epiproct, which was similar in both genders in the River Cova (Table 3). 


Table 3. Values of the comparison test of the variables between male and female exuviae of Boyeria irene 
in the Rivers Cova, Frio and Eresma rivers. HW: head width. BL: body length. PA: length of the paraproct. 
EP: length of epiproct. W+: greater width of the prementum. W-: smaller width of the prementum. PR: 
length of the prementum. F: ANOVA. H: Kruskal-Wallis test. P = probability. 


























Cova Frio Eresma 
HW | F=153.56 H = 31.026 F = 38.88 
P= 0.000 P=0.000 P=0.000 
BL | F=113.27 H=0.313 F = 4.76 
P= 0.000 P=0.576 P=0.033 
PA | F=18.61 H= 1.554 F = 4.76 
P= 0.000 P=0.209 P=0.033 
EP | F=2.79 F = 0.80 F = 3.47 
P= 0.097 P=0.374 P=0.068 
W+ | F= 102.37 F = 9.69 F= 16.14 
P= 0.000 P=0.002 P=0.000 
W- | F=70.12 F=1.01 F = 26.38 
P= 0.000 P=0.316 P=0.000 
PR | F=66.61 F = 10.025 F = 17.72 
P= 0.000 P=0.001 P=0.000 

















In males and females the values of the variables measured in the exuviae varied significantly between 
the three rivers (Table 4). 


Table 4. Comparison between the three rivers analysed in terms of the variables measured in exuviae of 
Boyeria irene. Each sex is analysed separately. H: Kruskal-Wallis test. F: ANOVA. P: probability. 





Male Female 
Test P Test P 

HW | F=32.78 0.00001 | H=48.60 | 0.0001 
BL | F=15.70 0.00001 | H=44.48 | 0.0001 
PA | F=7.26 0.001 H=18.02 | 0.0001 

EP | F=1441 0.00001 | H=19.10 | 0.0001 

GO H = 37.65 | 0.0001 
W+ | F=26.14 0.00001 | H=63.17 | 0.0001 
W- | F=17.21 0.00001 | H=55.17 | 0.0001 
PR | H=44.05 | 0.0001 H=81.58 | 0.0001 
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In the PCA analyzing the three rivers together, PC1 and PC2 component values for variance were: a) 
for River Cova 58.87% and 18.06% in males, 51.60% and 14.00% in females; b) for River Frio 56.62% 
and 20.86% in males, 59.74% and 14.63% in females; c) for River Eresma 57.70% and 23.35% in 
males, 56.90% and 15.31% in females. 


According to PCI values (Figure 1), epiproct is the variable with greater significance in males from 
the Rivers Cova and Eresma, while the length of the prementum had the greatest significance in 
males from the River Frío. In females, epiproct and prementum were the variables with greater sig- 
nificance in River Frío, and epiproct and paraproct in the River Cova. For the River Eresma epiproct 
and paraproct had similar values. The significance of the rest of the variables changed depending on 
the rivers. 
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Figure 1. Values of the loading of the variables in PCA 1 of male and female exuviae of Boyeria irene in 

each sampled river. HW: head width. BL: body length. PA: length of the paraproct. EP: length of epiproct. 


W+: greater width of the prementum. W-: smaller width of the prementum. PR: length of the prementum. 
AU variables were logarithmically transformed. 
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For males from the three rivers, the components PCI and PC2 (Figure 2) explained 59.07% and 
12.75% of the variance, respectively. With respect to PC1 all variables were located on the positive 
side, while with respect to PC2 the variables HW, and the three for the prementum (PR, W+ and W-) 
were on the positive side, and the variables BL, PA and EP were on the negative side. This implies 
that PC2 discriminates variables relating to the head rather than variables concerning the rest of the 
body. With respect to PC1 (Figure 2), on the negative side were found the exuviae from the River 
Frío, which have the lowest value for all the variables except BL. The exuviae from the Rivers Cova 
and Eresma were not different on this axis, because there was no clear difference between the two 
sites for the variables measured as a whole. Regarding PC2, the exuviae from the Rivers Cova and 
Eresma were clearly separate: on the positive side were found those from the River Cova, which have 
the largest values for HW, PR, W+ and W-, and the lowest values for BL and EP. 


0.20 


Component 2 





-0.20 
Component 1 


Figure 2. PCA of male exuviae of Boyeria irene in the Rivers Cova (grey points), Frío (red points) and 
Eresma (blue points). HW: head width. BL: body length. EP: length of epiproct. PA: length of the para- 
proct. W+: maximum width of the prementum. W-: minimum width of the prementum. PR: length of the 
prementum. All variables were logarithmically transformed. 


For females (Figure 3) the PC1 and PC2 components explained 59.49% and 16.53% of variance, re- 
spectively. With respect to PCI, all variables were located on the positive side. For PC2, the variables 
HW, PR, W+ and W- were on the positive side, as was also the case in the males, but they differed 
from the males in that BL was on the positive side. Finally, GO was also on the positive side in PC2. 
With respect to populations, the exuviae from the River Frio were again located on the negative side 
of PCI with the lowest values of all variables. The exuviae from the River Eresma were plotted 
around the intersection of the two axes, while those of the River Cova tended to the positive side on 
PCI and PC2, possibly influenced by the greater length of the prementum, which has a higher value 
in the River Cova than the River Eresma. 
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Component2 





Component 1 


Figure 3. PCA of female exuviae of Boyeria irene in the Rivers Cova (grey points], Frío (red points) and 
Eresma (blue points). HW: head width. BL: body length. EP: length of epiproct. PA: length of the para- 
proct. W+: maximum width of the prementum. W-: minimum width of the prementum 


Discussion 


Exuviae collected in the three sites fulfilled the previously established condition of differing in size 
among rivers. There were three variables (epiproct, paraproct and prementum) that stood out over 
others in the PCA and, therefore, could be considered representative in these populations. This was 
true for males and females. Furthermore, these three variables coincided in the three rivers analyzed. 
Therefore, in order to biometrically analyze B. irene populations in the Iberian Peninsula, the best 
variables analyzed here were the epiproct, paraproct and prementum. 


Epiproct, paraproct and prementum have been used previously in the taxonomies of odonates (see, 
among many others, Heidemann & Seidenbusch, 2002; Doucet, 2010; Tennessen, 2019). Their mor- 
phology and/or length aid differentiation between the species of various genera, above all in the way 
they join onto other body structures, such as in the genus Ladona (Tennessen et al., 2017). Often the 
length of the epiproct relative to the length of the paraproct is a great help for identifying the species 
correctly, as is the case with the genus 7ramea (Landwer & Sites, 2005) and other genera distributed 
in the Palearctic (Doucet, 2010). However, less attention has been paid to the intraspecific biometric 
variation in odonate larvae. In this sense, larval morphological characters vary with the habitat type 
occupied (Sahlén et al., 2008). Larvae living on the bottoms of streams, as is the case with B. irene 
(Velasco et al., 2018), may invest less in direct foraging organs (Giacomini & De Marco, 2008) and 
more in an effective system for reducing predation risk. In the populations analyzed in this study, the 
variables PR, EP and PA varied significantly between rivers, suggesting possible local adaptations to 
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environmental conditions. It is important to know the biometric differences between populations 
well, because variations in body size could play a role during competition for resources in adult and 
larval stages (Bowman et al., 2018; Minot et al., 2019). There is no doubt that knowing the morpho- 
logical differences between populations better would help us to improve our management and con- 
servation of the various odonates species (Tennessen, 2019). In the genus Boyeria, Müller et al. 
(2012) also refer to epiproct and paraproct as useful variables in interspecific analysis between B. 
irene and B. cretensis. Our data confirm that, in addition, these variables are valid for intraspecific 
analysis in this case, at least for B. irene, and that they could be used to analyze variations in popula- 
tions of this species. 


In the present study, the other variables analyzed in the exuviae of B. irene do not appear to be of 
value in this sense. Head width (HW) has been used to establish which instar B. irene larvae 
belong to (Ferreras-Romero, 1997), but the value of this variable changes on many occasions dur- 
ing ecdysis because the opening of the cuticle happens in an area of the body near the head (Cor- 
bet, 2004). This detracts from its value as a variable that is suitable for the biometric analysis of 
exuviae. Body length (BL) is an easy measurement to obtain for exuviae in situ, because it only 
requires a digital calliper. But there are other data that have significant importance in the set of 
variables when the three rivers are analyzed together. Furthermore, the intersegmental skins of 
the exuviae can be stretched out at emergence or not, and so the BL produces high variance. 
Consequently, BL is not suitable for the biometric characterization of B. irene populations, and 
should not be used in population analysis, at least for the Iberian Peninsula. The length of the 
gonapophyses (GO) is valid for separating Iberian and North African populations of Cordule- 
gaster boltonii (Casanueva et al., 2020), but in B. irene it is of little importance, since its value in 
the PCA is the smallest of all, and is therefore not useful for population analysis. The size of the 
gonapophyses in anisoptera larvae varies between species (Heidemann & Seidenbusch, 2002), 
but up till now this feature has not been used to study populations. On the other hand, adult fe- 
male genitalia of dragonflies have been widely used for taxonomies, phylogenetics, etc. (see, for 
example, Matushkina, 2011 and the references cited there). Clearly, more research is needed to 
establish whether the gonapophyses in larvae of different species is a valid structure for conduct- 
ing population research. Finally, the variables W-- and W- have been proposed for taxonomic 
study of larvae and exuviae in Anisoptera (Doucet, 2010, among others). However, Casanueva et 
al. (2017) and El Haissoufi et al. (2018) have shown that in C. boltonii these are not efficient vari- 
ables to use for intraspecific biometric analyses. Our results confirm that they should not be used 
to this end in B. irene either. 


Some authors have differentiated between genera and species of odonates by combining morpho- 
logical data and molecular analyses (Torres-Pachón et al., 2017; Baskinger et al., 2008; Hovmöller & 
Johansson, 2003). In B. irene Kohli et al. (2014) found that the populations present in NW Spain 
(Galicia) could be genetically differentiated from those of the rest of the Iberian Peninsula but these 
researchers did not compare biometric data. If the three measures that we propose in this work are 
valid for the population of Galicia, it is something that requires further research. 
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